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TECHNICAL MEMORANDUM X-G4939 


LIMIT CYCLE ANALYSIS OF LARGE SPACE TELESCOPE WITH 
CMG NONLINEARITY 

I. INTRODUCTION 


The purpose of this investigation is to study the existence of self-sustained oscilla- 
tions in the Large Space Telescope (LST) system due to the presence of nonlinear gimbal 
friction of the control moment gyroscope (CMG), if it is used as the momentum exchange 
controller for LST [1,2]. A continuous-data, single-axis model of the LST control sys- 
tem is considered with the CMG nonlinearity modeled as shown in Figure 1 . One of the 
major differences seen in the various dynamic models of the CMG is in the accurate 
representation of the nonlinear friction characteristics of the CMG. The running friction 
of the gimbal axis consists of tachometer brush friction and the hysteresis drag which is 
associated with the brushless dc permanent-magnet torque motor. These frictional torques 
do not have significant effect on large motion gimbal-rate control loop performance, other 
than the power loss required to overcome the drag torque. However, the accurate repre- 
sentation of these torques for small motion has a pronounced influence on the gimbal- 
rate control loop small-motion performance. Whereas these effects have been second order 
for most previously reported systems, their effect is significant for the highly accurate per- 
formance requirements of LST. Therefore, the accurate modeling of these torques is 
essential in predicting the existence of limit cycles and the control loop response to very 
low level commands. 

Using the solid friction model reported in Reference 3, a mathematical model of 
the CMG friction nonlinearity is derived for continuous describing function analysis. It is 
shown that the input-output characteristics of the nonlinearity are only a function of the 



Figure 1 . Single-axis model. 







amplitude A of the Input sinusoid and the nonlinearity parameters, A continuous describ- 
ing function N(A) is derived for the CMC nonlinearity, using the analytical frictional 
torque expression, and conditions for self-sustained oscillations are determined analytically. 
For the expected nominal values associated with CMG and LST parameter, an absence of 
limit cycle behavior is predicted. Digital computer simulations have been carried out to 
corroborate the analytical results. 


II. A SINGLE-AXIS MODEL OF THE LST CONTROL SYSTEM 


A simplified single-axis model of the LST fine-pointing control system is shown in 
Figure 1, The block diagram represents a simplified version of the LST control system 
obtained by neglecting the bending modes of the vehicle and the current loop and back 
emf of the CMG. However, the results of this report can be extended to more complex 
system models without difficulties. 

The transfer function G 5 (s) denotes the vehicle dynamics of the LST as a single 
inertia I. The dynamics of the CMG are also modeled by a single inertia with the trans- 
fer function G 4 (s). The transfer functions G 2 (s) and G ? (s) represent the CMG gimbal 
controller. The torsional feedback of the CMG output axis has been neglected. The block 
N denotes the CMG gimbal friction nonlinearity. The assumed characteristics of the non- 
linearity are described and mathematically treated in the next section. 

The controller for the LST vehicle is represented by the transfer function G, (s') 
which is of the proportional-plus-derivativc (attitude-attitude rate) type. Since the objec- 
tive of this report is not the design of the LST system, this conventional form of Gj(s) 
is chosen for the purpose of computer simulation. It can be shown that the controller 
Gj (s) is optimal in the sense of pole-placement design without the nonlinearity. Simula- 
tion results show that the system’s response is quite satisfactory with this controller, 

A list of symbols that describe the LST system variables and parameters, together 
with their numerical values, is given in Table 1. 


III. MATHEMATICAL MODEL OF THE CMG FRICTIONAL NONLINEARITY 


The running friction of the gimbal axis of the CMG consists of both the tachometer 
brush friction and the hysteresis drag, which is associated with the brushless dc permanent- 
magnet torque motor. These frictional torques do not have significant effect on large- 
motion gimbai-rate control loop performance other than through the power loss required 
to overcome the drag torque. In ordinary control system practice, these frictional charac- 
teristics can be either modeled as viscous or coulomb frictions or totally neglected if the 
amplitude of the frictional torque is relatively small. However, because of the fine-pointing 
accuracy and stability requirements of the LST system, accurate representation of the 
frictional torque becomes extremely important. For fine-pointing studies of the LST, the 
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TABLE 1. LST SYSTEM PARAMETERS 


Variable 

Description 

Magnitude 

H 

CMG Angular Momentum 

, 271 Nms 

I 

Vehicle Inertia 

1.354 X 10 s Nms 2 

*GE 

Gimbal Inertia 

5.0 Nms 2 


Giinbal Rate Loop Integral Gain 

1.354 X 10 4 Nm 


Gimbal Rate Loop Gain 

379 Nms 

Ao 

Vehicle Controller Coefficient 

2 X 10 4 (rad/s)' 1 

Ai 

Vehicle Controller Coefficient 

3 X 10 3 (rad)' 1 

e 

i 

Error input Command to CMG 

0.136 Nm 

! 

Tf 

Torque Output of Nonlinearity 

Nm 

T f 

*0 

Running Friction Torque 

0.136 Nm 

T 

Gimbal Bearing Parameter 

1.871 X 10 s (Nms)* 

5 

Gimbal Position 

rad 

5 

Gimbal Velocity 

rad/s 

0 

Vehicle Position 

rad 

6 

Vehicle Velocity 

rad/s 


accurate modeling or the frictional torques is essential in predicting the existence of limit 
cycles and the control loop stability when the system is subject to very low-level disturb- 
ances and commands, 

In this section, a mathematical model of the CMG friction nonlinearity is derived, 
based on the Dahl solid friction model [3,4] . The purpose of this model is that it can 
be used for the continuous-data describing function analysis of the nonlinear LST system. 








A mathematical and simulation model of solid friction, such as that encountered in 
a bearing, was reported in Reference 4. Subsequent experimental studies have indicated 
that the frictional characteristics of the CMG can be approximated by this solid friction 
model. A more accurate model that has the same general characteristics as the one used 
herein has recently been developed [51. However, since it is not precisely the same as 
that used in References 3 and 4, studies arc now underway at NASA’s Marshall Space 
Flight Center, by the co-authors of Reference 5, and other CMG developers, to develop 
an understanding of the effects of the more precise model on CMG and LST dynamics. 

If CMG output shaft friction torque versus shaft attitude (6) or attitude rate (t$) were 
plotted, it could be seen that the simplified model used herein yields the envelope of the 
plots obtained with the model of Reference 5. Hence, it is felt (but as yet unproven) 
that the limit cycle predictions yielded by using the simplified model are approximately 
correct. 


As indicated in References 3 and 4, the key to tiie simulation model is that the 
frictional torque, Tf, is a function of the CMG shaft displacement, 5, and can be differ- 
entiated with respect to time. Thus, 


dTf dT f dS • 
dt d5 dt f 


(I) 


The frictional torque Tf can be generated by integrating both sides of equation (1) with 
respect to t . 

in the simulation model, Tf is generated by a function generator with input Tf. 

The output of the function generator is multiplied by 5 to give the right-hand side of 
equation (1), which is then integrated to yield Tf. Figure 2 is a simulation diagram of 

the CMG friction nonlinearity. The input of the model is 6, and the output is Tf. The 
relays are necessary because Tf is sensitive to the sign of 5. 

it has been demonstrated experimentally that for solid rolling friction, the rela- 
tion between Tf and Tf may be approximated by a square-law expression, 

Tf 7(Tf-Tf Q ) 2 , T f < T ffl = 0 , T f > T fo I (2) 


where y is a positive constant and Tf o is the normalizing factor of Tf, However, 
the frictional torque is also velocity (5) dependent, as shown in Figure 2. Therefore, 
equation (2) should be written 



( 3 ) 
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Figure 2. Simulation of CMG friction nonlinearity, 

where 


T f . = T f sgn 5 . (4) 

Integrating both sides of equation (2), one obtains 

5 + C = -7 (Tf + - Tf ) 1 , 6 > 0 (5) 

and 

6 ■ C = -7.(T f “ - T^)* 1 , 5 < 0 j (6) 

where C is the constant of integration and 

T f = T f + , 5 > 0 

= T f ' , 6 < 0 (7) 


5 







If one assumes a sinusoidal input for 5 , such ns 


5 = A cos cot 


( 8 ) 


then 


C * A + [ 7 (R + I) T f J- 
R = -a' 1 + \J (a 2 + l)/a 2 


a = 27AT fo 


(9) 

(10) 

(ID 


and solving equations (5) and (6) for Tf , one obtains 


Tf = T f + = -t7(A cos cot + C)]*' + T fo 


-i 


5 > 0 


= Tf = -[7(A cos cot - C)r‘ - T f() , 8 < 0 


( 12 ) 


IV. A DESCRIBING FUNCTION FOR THE CMG GIMBAL 
FRICTION NONLINEARITY 


In order to derive a describing function for the assumed CMG nonlinearity, one 
assumes a sinusoidal input to the nonlinearity of the form of equation (8). The output 
of the nonlinearity, Tf , is approximated by the fundamental component of its Fourier 

scries. Since the input-output relation of the nonlinearity is symmetrical about the zero- 
torque axis, there is no dc component in the scries. Thus, 


Tf « T|(t) = A| sin cut + Bj cos cut , 


(13) 


where 


2lT 


A, « 


f Tf sin cot dcot = - f Tf sin cot dcot + 


2tt 

J Tf + sin cot dcot 

7 T 


( 14 ) 
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and 


B, 



2tt 

J Tf cos cot dcot 
0 


l 

7T 


2 rr 

f Tj-' cos cot dcot 
0 


+ 


i 

7T 


27T 

J* Tf + cos cot dcot 


( 15 ) 


Equation (13) can also be written as 

T, (t) = Va, 2 + B, 2 cos (cot - 0) 


06 ) 


where 


0 = tan’ 1 (Aj/Bi) 


(17) 


or in phasor notation, 

T i (ico) = >/ Aj 2 + B| 2 exp [i tan" 1 (Ai/B,)] . (18) 

The describing function is given by 

N(A) = T, (ico)/5(ico) = (B,/A) - i(A,/A) = N, + iN 2 . (19) 


The coefficients A, and B( can be obtained by substituting equation (12) into 
equations (14) and (15), respectively, and performing the appropriate integrations. The 
results are: 


A, = -(4T fo /tt) + 2Cn[(C + A)/(C- A)1/7 tA T 

= (2T f M[fin(a + y/ a 5 + l) 2 - 2a] (20) 
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and 


B, = 2[(CA/ C 2 - A 2 ) - i \lyA 

■ (4T f /o) >/[a 2 +a + I + (a + \)y/ a* + l]/2(a +y/a 1 + 1) - 1 . (21) 


For stability analysis, it is of interest to determine the behavior of *I/N(A). This 
may be computed readily by use of equations (19) through (21), Figure 3 shows the 
magnitude (db) versus phase (degrees) plots of -I /N(A) for y- 1.871 X 10 s , 1.871 X 
1 0 s , and 1.87 X 10 7 (Nms)* 1 , as the magnitude of A varies. 



Figure 3. Plots of -i/N and G(S). 


Figure 3 shows that as A aproaches zero, the magnitude of -1/N(A) in db 
approaches 20 logi 0 [ l lyV f Q ] and the phase is -180 degrees. This asymptotic behavior 

of -1/N(A) for small values of A may be derived analytically be expanding the 
logarithmic term in equation (20) and using relations (7) through (9): 
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Hm (Aj/A) = 0 
A->0 


( 22 ) 


Similarly, by expanding the square root term of equation (21) into a power series one 
obtains 

lim (Bi/A) =o , (23) 

A->0 

where 

= ? T ? 0 

Thus, 

lim [-1/N(A)1 = lim {-[(B,/A) - i(A|/A)]}'‘ = -1/a 
A->0 A->0 

For large values of A , the value of C becomes 

lim C = lim (A + 2/yT^) = lim [A] -+ 00 
A-x» A-**> A-** 3 

Then, 

lim (A, /A) = lim (-4Tf o /7rA) - -0 , (27) 

A-x» A-*» 

and 

lim (B, /A) = +0 (28) 

A-Hpo 


(24) 


(25) 


(26) 
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lim t-l/N(A)l « -i/0 = oo/-270° 

A-K50 


( 29 ) 


As shown in Figure 3, the gain-phase plots of -l/N(A) approach °°/-27 8° as A-+ f >° 
lor all values of y and T|- q . 

Alternately, a general map of N 2 versus Nj (as a function of the argument A) 
may be developed for use in this and future limit cycle investigations for this particular 
form of nonlinearity (Fig. 4) [6] . It is universal in the sense that it permits one to choose 
any design values for Tf and 7, and sec their effect on the N 2 versus Ni map. 

N, 


a 



NOTE; describing FUNCTION ASYMPTOTES: — 
VALUES FOR A ARE FOUND FOR CURE 
CURVES No. t It No. 2, RESPECTIVELY, 

BY READING THE VALUE OF Nj/o FOR THE 
POINT IN QUESTION AND THEN APPLYING 
THE APPROXIMATION INDICATED. 

Figure 4, Normalized Dahl describing function locus. 


V. STABILITY ANALYSIS BY USE OF THE DESCRIBING FUNCTION 


The system described in Figure 1 appears to be suitable for describing function 
analysis because it is low-pass and the system parameters are assumed time-invariant. The 
condition for a self-sustained oscillation in a linear system with a nonlinearity is given by 


1 


1 + N(A) G(ico) = 0 


(30) 




where N(A) is the describing function of the nonlinearity and G(iw) is tho transfer 
function of the system which is seen by the nonlinearity. The solution of equation (30) 
can be obtained graphically by plotting the G(ico) curve for the system and the -1/N(A) 
curve for the nonlinearity. The point of intersection of these curves will yield the solu- 
tion to equation (30). The magnitude and frequency of the oscillation corresponds to 
the values of A and co at the solution point. 

For the LST system in Figure 1 , the transfer function which the nonlinear element 

sees is 


G(s) 


Is 3 , 

llcgs 4 + IKpS 3 + IKjs 2 + KjHAiS + K|HA 0 


(31) 


The relevant portions of the frequency-domain plot of G(s) are super-imposed 
on the -1/N curves of Figure 3 in db versus phase coordinates. With y = 1.87 X 10 7 , 
the -1/N curve intersects the G(s) curve at two points. The stable point for sustained 
oscillations is the one on the left at the higher frequency. The approximate magnitudes 
for frequencies of the oscillations arc 3 X KT 7 rad and 1.8 rad/s, respectively, for the 
unstable limit cycle, and 2 X 10' 5 rad and 5.6 rad/s, respectively, for the stable limit 
cycle. The curves in Figure 3 also show that for y considerably smaller than 1.871 X 
10 7 , the system will exhibit a stable response. 

Following the alternate approach begun in the previous section, limit cycle condi- 
tions and characteristics may be predicted, using the parameter plane technique [6] . 
Specifically, one may place the describing function locus (Fig. 4) on an Ni , N 2 param- 
eter plane stability map to analyze and predict limit cycle behavior. The result is por- 
trayed on Figure 5 and agrees with the stability analysis technique portrayed on Figure 3. 


VI. SIMULATION 


To corroborate these results, the LST system was simulated on a digital computer. The 
CMG friction nonlinearity was simulated as shown in Figure 2. For the computer simula- 
tion, the input to the LS’i system, T c , was set to zero, along with all of the initial stater , 
except for the vehicle position, Q , which was set at 5 X 10‘ s rad. This value was chosen 
so that the input signal to the nonlinearity, 8 , would be large enough to cause the torque 
to saturate but not so large as to exceed the limiting value of the input signal. The 
quantities plotted from the simulation runs are vehicle attitude and attitude rate, 6 and 
0, and friction torque. Tf. 
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Figure 5. Nj - N 2 parameter plane. 


Figure 6 shows the plots with 7 = 1.871 X 10 7 . It is observed from the plot of 
Tf in Figure 6 that the system has a sustained oscillation. This oscillation is not seen on 

the other plots because of the large initial transients. For this reason, the vertical scale 
of the plot is adjusted (Fig. 7) to show the oscillations. The frequencies and magnitudes 
of oscillations obtained are quite close to the predicted values. The small discrepancy is 
partly attributed to the quantization caused by the nonlinearity implementation on the 
digital computer. 

Figure 8 shows the response plots with 7 = 1.871 X 10 s . As predicted, the 
system is stable for the lower values of 7 . 
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TIME (t) 


Figured. Simulation (y = 1.871 X 10 7 ). 


VII. CONCLUSION 


For numerical values considered to be representative of the LST and CMG’s (if 
they are used), analysis of the simplified continuous-data model with a solid friction CMG 
nonlinear characteristic indicates (and analogue simulation confirms) the absence of limit 
cycle behavior due to the CMG nonlinearity. It must be borne in mind that the analysis 
has not yet included the dynamic effects induced by the sampling phenomena of the 
onboard digital computer. Experience with past spacecraft would lead one to expect 
quantization, for example, to induce limit cycling of the plant. 
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Tf{X10 2 )(Nm) 0(X1O‘ s )(rad/s) 0(X1O' 7 )(rad) 


0.04 



TIME (i) 



TIME (i) 



Figure 7. Simulation (7 = 1.87! X 10 7 ): steady state. 
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